Electrically driven magnetization of diluted magnetic semiconductors actuated by 

Overhauser effect 
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It is well-known that the Curie temperature, and hence the magnetization, in diluted magnetic 
semiconductor (DMS) like Gai-^Mn^As can be controlled by changing the equilibrium density of 
holes in the material. Here, we propose that even with a constant hole density, large changes in 
the magnetization can be obtained with a relatively small imbalance in the quasi-Fcrmi levels for 
up-spin and down-spin electrons. We show, by coupling mean field theory of diluted magnetic 
semiconductor ferromagnetism with master equations governing the Mn spin-dynamics, that a mere 
splitting of the up-spin and down-spin quasi-Fermi levels by O.lmeV will produce the effect of an 
external magnetic field as large as IT as long as the alternative relaxation paths for Mn spins (i.e. 
spin-lattice relaxation) can be neglected. The physics is similar to the classic Overhauser effect, also 
called the dynamic nuclear polarization, with the Mn impurities playing the role of the nucleus. We 
propose that a lateral spin-valve structure in anti-parallel configuration with a DMS as the channel 
can be used to demonstrate this effect as quasi-Fermi level splitting of such magnitude, inside the 
channel of similar systems, have already been experimentally demonstrated to produce polarization 
of paramagnetic impurity spins. 

PACS numbers: 



I. INTRODUCTION 



Electrically driven magnetization of diluted magnetic 
semiconductors (DMS) has the potentiality to open up 
new avenues on the map of magneto-electronics and spin- 
tronics [H, • I n this regard, electrical manipulation of 
magnetization has already been demonstrated an d 
theoretically proposed [l(J [ll[ . The Curie temperature, 
in these methods, were controlled by changing the car- 
rier concentration (Fermi level) while keeping the car- 
rier spin-subsystems in an equilibrium among themselves, 
that is, keeping the quasi-Fermi levels for up-spin (/.*-]-) 
and down-spin (/ij.) carriers equal. In contrast, in this 
work we propose that even with a constant carrier den- 
sity, large changes in the magnetization can be obtained 
with a relatively small imbalance in the spin population, 
that is, a small difference in /i-j- and fi^. We also propose 
a structure (fig.[lja)) for demonstrating the effect that is 
within current experimental capabilities. In essence, our 
proposed scheme is similar to the optical manipulation 
of magnetization in refs. fl2l ] and where an imbalance 
in the spin-population is attained by shining circularly 
polarized light. 

Our proposed effect represents a non- equilibrium mag- 
netization resulting from a non-equilibrium bath (the car- 
rier spins) constantly trying to restore equilibrium via 
spin-flip process due to exchange interaction with local- 
ized spins which gets polarized in the process. Indeed 
the physics is similar to the classic Overhauser effect, 
also called the dynamic nuclear polarization (DNP) fl3l ). 
with the Mn impurities playing the role of the nucleus. 



To our knowledge this effect has not been employed to 
actuate non-equilibrium magnetization by electrical exci- 
tation although magnetization by optical excitation via, 
possibly, the same effect 0, H3 and demagnetization via 
the opposite effect [13] has been experimentally observed. 
Our proposed effect involves electrically driven dynami- 
cal polarization of interacting spins (where the polariza- 
tion of a particular localized spin is affected by the po- 
larization of the neighboring localized spins) and, hence, 
would be an extension of similar effect studied in the 
context of non-interacting spins f]~5l— Tl7| . Due to this ef- 
fect a splitting of the up-spin and down-spin quasi-Fermi 
levels in the channel (fig. [Ua)) by O.lmeV can have the 
same effect as an external magnetic field of IT (fig.QJb)). 
Splitting of this order can be attained by spin-injection 
into semiconductors [HI and has recently been demon- 
strated in an n-channel GaAs lateral spin-valve operated 
with anti-parallel contacts to actuate dynamical polariza- 
tion of non-interacting spin [l!| [2(| ■ A similar p-channel 
Gai-^Mn^As (x ~ 0.05) structure should be suitable for 
the demonstration of the proposed effect. 



II. MODEL OVERVIEW 



'Electronic address: lsiddiquOpurdue.edu 



A number of theoretical papers [2lM27| have modeled 
the appearance of ferromagnetic ordering among the Mn 
ions in Gai-^Mn^As, interacting via the itinerant holes, 
in terms of a mean field description and explains the 
experimentally observed pel [29| temperature variation 
of magnetization in these materials. We adopt exactly 
the same model as ref. [2(| and have modified it to take 
into account the non-equilibrium aspect by: i. intro- 
ducing two different quasi-Fcrmi levels for up-spin and 
down-spin holes (/i-f and ^tj, respctively) inside channel, 
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and ii. writing a master equation to describe the non- 
equilibrium dynamics of the Mn spins that was used in 
ref. [2(i| to semi-quantitatively explain the experimental 
observation of Mn spin-dynamics in ref. I19i . In essence, 
our model is the same as that in refs. 3 3(|, which also 
studies the magnetization dynamics of DMS and uses a 
more sophisticated valence band description. 

With /j-j- — /.ii = we get essentially the same results as 
ref. [26[ (solid curves in figs. [TJb) andQJc)). Under non- 
equilibrium situation (fj^ — ^ 0), according to our 
model, we expect to see a strong ferromagnetic ordering 
among the Mn ions (fig.[l]) for moderate values of fjrf —f^i 
due to reasons that we will discuss later in this paper. 
A moderate magnitude of /i-f — ~ 0.1 meV is quite 
feasible inside the channel of an anti-parallel lateral spin- 
valve structure (l9l . |20| and can be understood in terms 
of a circuit model presented in ref. jl9j to explain the 
experiment therein. We use the same model later in the 
paper to estimate /if — /ij.. 



III. THEORY 
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FIG. 1: (color online) . (a) Schematic structure: A dilute mag- 
netic semiconductor (DMS) Gai-zMn^As connected through 
nonmagnetic (NM) spacers to two ferromagnetic (FM) con- 
tact in antiparallel (AP) spin-valve configuration. Under bias 
V , an electronic current /, injecting spin-polarized carriers 
(holes), flows, (b) Mn magnetization M, scaled by Mn satu- 
ration magnetization Mo, vs. temperature T, scaled by the 
Curie temperature T c , of the channel in (a) for different quasi- 
Fermi level splitting /i^ — /tj. (Mt(4.) being the channel quasi- 
Fermi level for up(down)-spin carriers) with external magnetic 
field B ex t — 0; inset: under different B ex t at equilibrium 
(/if — fii—O). (c) Hole magnetization m, scaled by satura- 
tion hole magnetization mo, vs. temperature T, scaled by the 
Curie temperature T c , of the channel in fig. QJa): under the 
same conditions as in fig. 0Jb); inset under the same condi- 
tions as in inset of fig. UJb). The parameter values for these 
calculations are: nMn ~ 5.0 x 10 20 cm~ 3 , n^/nA/n = 0.08, 
m* = 0.5m e , ao — 5.65 A, and J — leV. 



The spontaneous ferromagnetic ordering of the local- 
ized Mn spins in a DMS material arises due to the hole 
mediated exchange interaction between them. In the con- 
text of mean field theory of DMS [26[ , the carriers 'feel' 
an exchange field due to the polarized Mn spins in addi- 
tion to any external magnetic field B ext , which separate 
the up-spin band from the down-spin band (fig. Ufa)) in 
energy by, 



(ex) 



= Jaln Mn (Sf n ) 



(1) 



where, ao is the lattice constant, gh is the g-factor of the 
carrier (hole), hb is the Bohr magneton, A( ea .) is the sep- 
aration between up-spin band and down-spin band due 
to exchange field and (S^) is the average z-component 
of S = 5/2 Mn spins: 



(2) 



F s being the probability of a Mn spin being in S^ In = s 
state (s = 5/2, 3/2, . . ., -5/2). 

The carrier band splitting in addition to splitting of 
quasi- Fermi levels for different carrier spins (fig. [DJa)) 
lead to a non-zero average z-component of hole-spin (S^) 
due to unequal carrier concentrations for up-spin (ji/i,f) 
and down-spin (n^.i), which arc given by: 



(si) = 77 



1 n/, )T - nh.x 



2 n ft ,t + n h . x 

dE(l - f m) (E))D m (E) 



(3) 
(4) 



where, f-^-^(E) is Fermi function with Fermi 
level fJ-f(i) an d temperature T: ff(i) (-^0 = 
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(a) 



FIG. 2: (a) Schematic carrier (hole) bands of the channel 
in fig. QJa) for different spins: E v ^^, /U-|-(4.) and A being the 
valence band edges for up(down)-spin holes, their quasi-Fermi 
levels and their band splitting respectively, (b) Energy levels 
and relaxation of Mn spins inside the channel in fig. Ola): 
Sz In denotes the six spin states of the S = 5/2 Mn spins, 
e denotes their energy difference, r-j-m-j-) denotes the rate of 
transition from lower (higher) z-component state to a higher 
(lower) z-component state and, n = —5/2, —3/2, . . ., 3/2. 



tcrnal magnetic field, which lead to energy level splitting 
(fig. lib)), 



6 — £(ei) + 9MnHBB ex t 

£(ex) = Jal n h(Sz) 



(5) 



(a new variable et ex ) appears in the above equation) 
where, gMn is the g-factor of Mn spins and £( ex ) is the 
splitting of Mn spin levels due to exchange field. 

Up to this point, all the equations and quantities, ex- 
cept different quasi-Fermi levels fjrf and /xj,, are essentially 
the same as ref. [26[ and takes into account the non- 
equilibrium effects of Mn spin-polarization on the holes. 
To consider the non-equilibrium effects of the holes on 
the Mn spins and the resulting non-equilibrium dynam- 
ics we solve the dynamic rate equation, which was also 
used in ref. [2(| Ho| , for the occupation probabilities F s : 



[l + cxp {(E - fi ni) ) /k B T}] , and £> t m is the 
three dimensional density of states for up(down)-spin 
carriers calculated assuming a parabolic band with 
an effective mass of m* and a band-edge at E v ^^ 
(fig-EJa)). E v -f, E v i and A are related by E v ^ — E v ^ = A 
and the charge neutrality condition: = rihrf + w>h,x- 

The Mn spins on the other hand feel an exchange field 
due to the spin-polarized carriers in addition to any ex- 
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at steady state (by setting dF s /dt = 0) under the nor- 
malization constraint: ^ s F s = 1 where, 
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(7) 



Tfi and (fig. E£b)), as valence band holes surround- 
ing the Mn spins act as their spin-bath and also as Mn 
spins have negligible spin-lattice relaxation rate in com- 
parison [l4j . are given by: 

T n = ^ J 2 J dED x {E + e)(l - h{E + e))D t (E)f t (E) 

(8a) 

T it = HE J 2 J dED t (E)(l - ff(E))Di(E + e)f±(E + e) 

(8b) 



IV. RESULTS: NONEQUILIBRIUM 
MAGNETIZATION 

The eqs. El El EJ El and Q] are solved sequentially and 
self-consistently (fig. E|) to calculate the Mn spin po- 
larization and hole polarization in figs. [TJb) and [TJc) 
respectively. Magnetization of the Mn and holes in- 
side channel material (Gax^Mn^As) is then calculated 
from: M = 2M (S* In )/5 and m = 2m (S%) respec- 
tively, where the corresponding saturation magnetiza- 
tions, M = 5gMnMsnMn/2 and m = g h [i B n h j2. The 
iterative loop in fig. [3] embodies a positive feedback loop 
that gives rise to equilibrium magnetization below Curie 
temperature in such materials. Such positive feedback in 
combination with Overhauser effect (fig. [3]) is what gives 
rise to non-equilibrium magnetization even above Curie 
temperature. The results of the calculations, shown in 
figs- Hlb) and HJc), use realistic parameter values [26| 
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and gives equilibrium magnetization characteristics (cal- 
culated by setting fi-t — = 0) similar to experimental 
observations [HI, H^j ■ Two sets of calculations were per- 
formed: 1) calculations for different values of ^ — /xj, by 
setting B ext = 0, whose results are shown in the main 
plots of figs.QJb) and [ITc), and 2) calculations for differ- 
ent values of B ext by setting fj^ — = (equilibrium) , 
whose results are shown in the inset plots of figs. HJb) 
and QIc). By comparing the main plots in figs. [TJb) 
and rjjc) with their corresponding inset plots, we ob- 
serve that the values of fj^ — /Uj, corresponding to the 
curves in the main plots maintain a proportionality rela- 
tion with the values of B ext corresponding to the similar 
curves in the corresponding inset plots (O.lmcV: l.OmeV: 
10.0meV= IT: 10T: 100T). Such proportionality rela- 
tion is, by no means, a coincidence and is maintained 
in the calculations done with different values of J, riMm 
and rih (results not shown in this paper). Although the 
Curie temperature, the saturation magnetization and the 
shapes of the magnetization vs. temperature curves are 
different for different values of J, riMni and n h due to 
the dependence of the exchange field on these param- 
eters (eqs. Q] and [5]) the strength of the effect remains 
the same, i.e. the magnetizations for a quasi-Fcrmi level 
splitting of O.lmeV without any external magnetic field, 
when alternative Mn spin relaxation paths can be ne- 
glected, is equal to the magnetizations for an external 
magnetic field of IT at equilibrium. As a result, chang- 
ing J, riMm and rih by changing the doping, changing the 
Mn mole fraction, introducing disorder [3l| or by using a 
different DMS material (that has itinerant carrier medi- 
ated exchange interaction of localized spins) will not play 
any significant role as far as the strength of the effect is 
concerned for reasons to be discussed in the next section. 



V. DISCUSSION 

The results in figs.QTh) and [1(c) can be anticipated if 
the effect of ^t-j- — /ij, is considered as an effective external 
magnetic field as far as the hole spin-polarization and Mn 
spin-polarization arc concerned. We can show that the 
functional dependence of (S^) on A and /i^ — /ij, obeys 
the following relation (see appendix lAf : 

(S, h )(A,Mt-/H) = (^(A + Mt-^,0) (9) 

The first term is the hole spin-polarization for a valence 
band splitting of A and a non-zero quasi- Fermi level split- 
ting fj^—fj,^ while the second term corresponds to the hole 
spin-polarization for an additional valence band splitting 
of /i^ — /i| over the orginal value A and a zero quasi- 
Fermi level splitting (equilibrium). At this point, one 
can observe that the results in the main plot (i.e. for 
Bext = 0) and the inset plot (i.e. for ^ — = 0) 
of fig- Die) (i- c - hole spin-polarization) correspond to 
(^(A^Mt-/^) i.e. (^)(A (ex) + Mt -/ir,0) (eq.ij) 
and (St? )(A( ea ,) + ghl^BB ex t,0) respectively. The mathe- 
matical equivalence of the last two expressions suggests 
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FIG. 3: Nonequilibrium magnetization arises from an inter- 
play of dynamical spin-polarization (Overhauser effect) that 
is driven by the quasi-Fermi level splitting — ^tj. between 
up-spin and down-spin holes and a positive feedback that is 
responsible for equilibrium magnetization below Curie tem- 
perature 



that, as far as hole spin-polarization is concerned, ^ — 
has the equivalent effect of an effective external magnetic 
field of B^ xt = (fif — fii)/(ghfJ.B)- On the other hand, for 
Mn spins it can also be shown that the functional de- 
pendence of (S^ In ) on e and fif — obeys the following 
relation (see appendix iBj) : 



(5f")(e, Mt -^) = (5f n )(e + Mt-/X4,0) (10) 

In this case, the first term is the Mn spin-polarization for 
an energy level splitting of e and a non-zero quasi-Fermi 
level splitting of ^ — /ij, while the second term corre- 
sponds to the Mn spin-polarization for an additional en- 
ergy level splitting of /i-j- — /i^ over the orginal value e and 
a zero quasi-Fermi level splitting (equilibrium). At this 
point, one can observe that the results in the main plot 
(i.e. for B ext = 0) and the inset plot (i.e. for ^ — ^ = 0) 
of fig- IUb) (i.e. Mn spin-polarization) correspond to 
(5f")(e (ea) ,Mt-|"i)i-e. (Sf n ){e {ex) +^- H , 0) (eq.HDJ 
and {S* ln )(e( ex ) + g h n B B ext ,0) respectively. The math- 
ematical equivalence of the last two expressions shows 
that, as far as Mn spin-polarization is concerned, ^ — /i^ 
has the equivalent effect of an effective external magnetic 
field of = {^—^i) I {gMn^B) acting on the Mn spins 
at equilibrium. For g^ ~ gi\i n ~ 2 the equivalent exter- 
nal magnetic fields for hole spins and Mn spins (B^ xt and 
B^ x ™), that are mentioned above, would be eqaul and be 
given by 



B 



eff _ Mt ~ Mi 



ext 



(11) 



Herein lies the strength of the effect: a mere difference of 
O.lmeV between /x-j- and ^tj, is strong enough to produce 
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the effect corresponding to that of an external magnetic 
field as large as IT. One can notice that the arguments 
presented above do not depend on J, tiMn or n h and 
require that the alternative relaxation paths for the Mn 
spins can be neglected (which entered through the ne- 
glect of spin-lattice relaxation rate while writing down 
the eqs. [8] and was subsequently used in the derivation of 
eg. 1 1 0[) . As a result, the strength of the effect ( eg. \11} ) is 
insensitive to a change in J, nMn, and nh, and, hence, 
is relatively insensitive to our choice of mean field ( eqs. Q] 
and\3j) as long as the alternative relaxation paths for the 
Mn spins can be neglected. Since 1) our central result 
(eg. Hip , originating from non-equilibrium magnetization 
dynamics, is relatively insensitive to our choice of mean 
field (eqs. Q]and E|), 2) eqs. [T] and [5] describe the equilib- 
rium temperature dependence of magnetization, at least, 
qualitatively (2(|, and 3) a model [3(|, which essentially 
uses our choice of mean field, has been used to success- 
fully explain recent experiments (l4j on non-equilibrium 
magnetization dynamics we feel justified in leaving it to 
future work to assess the need for improving eqs. Q] and [5] 

Experimental realization of quasi-Fermi level splitting, 
similar to the values mentioned above, has already been 
demonstrated to drive dynamical polarization of non- 
interacting spins [ill, H(| with a structure similar to the 
one shown in fig. [T{a). One can quite legitimately en- 
vision that with further improvement of spin-injection 
process in terms of contact polarization P c and contact 
conductance and hence the parallel terminal conductance 
G|| one can achieve even higher ^ — (i^ leading to higher 
degree of ferromagnetic ordering that would otherwise 
require immensely large magnetic field (fig. [T]). The ma- 
terial property of the DMS that acts against attaining 
quasi-fermi level splitting between up-spin and down-spin 
holes is the spin lifetime of valence band holes t so in- 
dependent of Mn that give rise to the spin-flip conduc- 
tance g so . The effect of all these ingredients of a spin- 
valve structure on ^ — (i^ can be concisely pictured in 
terms of the circuit model in fig. HJa) which we have 
adopted from ref. [l9| and is valid for a channel length 
in fig. [l]Ja) that is smaller than the spin-diffusion length. 
Since the spin-diffusion length of the magnetic semicon- 
ductors has not been reported in the literature to the best 
of our knowledge we arc unable to conclusively comment, 
as far as the channel length is concerned, on the scope 
of the analysis to follow and will limit our analysis to 
the thinnest possible (2D) channel having a thickness of 
atomic dimensions in the transport direction. Neverthe- 
less, it will touch upon some key ingredients that affect 
/i-f — Moreover, the derivation of eqs. [5] and [TU] do not 
rely on the shape of the density of states as long as the 
density of states for both up-spin and down-spin carriers 
have the same energy dependence (so that one can write: 
D-f-(iJ) = Di(E + A)). As a result, as far as the strength 
of the effect is concerned, whether the channel is 2D or 
3D does not play any significant role. 

Upon simplification of that circuit model as shown in 
fig. Hb) we get: ^ - H = qVP c {\ + f^)" 1 where, V 
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FIG. 4: (a) Circuit model of the strcuture in fig.QJa) (adopted 
from ref. [3): g so , G?t(4.).i) Gt(j.),2, 9a and, gp are the car- 
rier spin-flip conductance independent of Mn spins, contact 
1 conductance for up(down)-spin, contact 2 conductance for 
up(down)-spin, majority spin conductance and, minority spin 
conductance respectively, (b) Simplified circuit diagram of 
the model circuit in (a): G\\ and P c = {g a — gp)/(ga +9fs) are 
the terminal conductance for parallel (P) configuration and 
the contact polarization respectively. 



is the applied bias and q is the electronic charge. We 
estimate g so ~ 10 10 f2 _1 m~ 2 for Gai-^Mn^As from the 
relation g so — ^-D^- using t so = lps (spin life-time 
of GaAs valence band holes (HJ) and a 2D density of 
states value D <~ 10 37 J~ 1 m -2 (estimated using the va- 
lence band effective mass used in fig. Lljb), which is of 
the same order of magnitude as the number calculated 
from 3D density of states for a thickness of atomic di- 
mensions (~ lnm)). For such value of g so if we use the 
terminal conductance and contact polarization values of 
refs. [H,[20| {G\\ - lf/fTW 2 , P c - 0.5) we estimate 
Mt ~ AH ~ lmeV for an applied voltage of V = IV. How- 
ever, for a t so = 10fs (the value used in ref. [lH) one 
would have to increase G|| to - lC^frW 2 (similar to 
the values in ref. [33j]) to get the same effect. 



VI. SUMMARY 

In summary, we have proposed a novel mechanism of 
attaining and controlling ferromagnetic ordering in DMS 
materials. We argue that the non-equilibrium accumula- 
tion of carrier spins drives the spin alignment of magnetic 
impurities in DMS via the exchange interaction which is 
responsible for the origin of the ferromagnetism in such 
materials in the first place. We find the effect to be 
quite strong when we consider that such degree of fer- 
romagnetic ordering would otherwise have to be attained 
by applying a very large magnetic field and, also that 
this strength is insensitive to a change in several crucial 
parameters (magnetically active Mn concentration, hole 
concentration, and strength of exchange interaction be- 
tween Mn spins and valence band holes) that describe 
the equilibrium ferromagnetism of such systems. We 
believe that, with the existing experimental sophistica- 
tion achieved over the years, our proposed scheme would 
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be realized in the near future and will usher in a new 
paradigm in the experimental investigation of ferromag- 
netic phase transition and also in the applications bene- 
fitting from the strong control of magnetism such as the 
magnetocaloric applications }34| . 

This work is supported by the Office of Naval Research 
under Grant No. N00014-06-1-0025. 



Appendix A: Derivation of eq. [9] 



Starting from eqs. [3] and |4] we get 

1 



(^>(A, Mt - Mi ) 



2n h 



J dE{D^E)-D x {E)} 



J dE{fr(E)D t (E) - f±(E)D±(E)} 



(Al) 



where was assumed to remain constant (varying with 
neither A nor fj,^ — /xj,) on the ground of charge neutrality. 
Now, for a given valence band splitting A and a given 
quasi-Fermi level splitting <fyt = fif — fj,^ (fig. [2]) we can 
write 

D t (E) = D l (E + A)=D(E) (A2) 
f t (E) = f t (E - Sfi) ee f(E) (A3) 

Substituting the above results into eq. IA1I 

(SH)(A, K-to) = J-\[dE {D(E) D(E A)} 
2n h [J 

J dE {f(E)D(E) - f(E + 5fi)D(E - A)} (A4) 
The part 

J dE {D(E) — D(E — A)} 
= J dED(E) - J dED(E - A) 

= J dED(E) - J dE'D(E') 

by performing a change of variable on the second inte- 
gration in the previous step. It, finally, leads to 

J dE {D(E) - D(E - A)} = 



Substituting the above result in eq. IA4I we get 

1 



(^}(A, Mt -^) 



2n h 



dE{f(E)D(E) 



- f(E + 8n)D(E-A)} (A5) 

Proceeding further, 
(^}(A, Mt - m ) 

= -L dEf(E)D(E) - J dEf(E + Sn)D(E - A) 
= 2^ {/ dEf(E)D(E) - J dE'f(E')D(E' - A - 8 ft) 



by performing a change of variable on the second inte- 
gration in the previous step. Finally, it leads to 

(^}(A, Mt - Mi ) = JL J \dEf{E)D{E) 

- f(E)D{E-A-5 l i)} (A6) 

while substituting //-(• — fj,± ee 5/j, = and A = A' + Sfi in 
eq. IA5I we find 



(SjXA' + ty.O) - ±-J d E{f{E)D{E) 

- f(E)D(E - A' - Sfi)} (A7) 

which, trivially, leads to 

(S£)(A + Mt-^,0) = -LJdE{f(E)D(E) 



- f(E)D(E - A - S/i)} (A8) 



The right hand sides of eqs. IA6I and IA8I being equal, we 
equate their left hand sides and arrive at the results in 
eq.El 



Appendix B: Derivation of eq. 1101 



Starting from eqs. |5] and [3 at steady state, we get 



F_ 



5/2 



F_ 



3/2 



F 4 



3/2 



F_ 



3/2 



F_ 



1/2 



'+5/2 



where, we have defined 



= iil 



(Bl) 



(B2) 



From eqs 
we get 



and probability conservation: F s = I 



5/2 



F_ 



3/2 



n" 



+5/2 



1 - a 

I - a 6 



(B3) 



Substituting the above results in eq. [2] we find 



(5 



Mn\ 



2(1 - a 6 ) 



(5a 4 + 8a 3 + 9a 2 + 8a + 5) (B4) 



which shows that (S^ In ) at steady state is solely a func- 
tion of a. Now, from eqs. 151 lA2l lAll and lB2l 

_ [ dED(E) {1 - f(E)} D(E + e - A)/(£ + e + fy) 
" ~ / d£Z?(£ + e - A) {1 - f(E + e + 6^)} D{E)f{E) 

(B5) 

For a given value of E the ratio of the integrands in the 
above equation 

D(E) {1 - f(E)} D(E + e - A)/(£ + e + gg) 
D{E + e- A) {1 - /(£ + e + <5//)} D{E)f{E) 

l-f{E)\ f l-fjE + e + S^ y 1 



= exp 



/(£) J I f{E + e + 5n) 
e + 5fi ' 
k B T 



7 



Using the above result while treating the integrations in From the above relation it trivially follows that 

cq. IB5I as summation over energy and making use of the 

identity 

Ni _ N 2 _ _ Ni _ _ T,i N i a ( e >MT-W-) = a '( e + M- -AH,0) (B6) 

From the above relation and the relation in eq. IB4I (which 
shows that (S^ In ) at steady state is solely a function of 
a) we arrive at the result in eq. [10] 



we finally get 



exp 



k B T 
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